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Abstract. We present high angular resolution radio snap-shot observations of seven nearby low-luminosity active 
galaxies (LLAGN) from the NUclei of GAIaxies (NUGA) survey. The observations were conducted with MERLIN 
and EVN/VLBI at 18 cm and 6 cm. At all observed angular resolutions and frequencies, we find indications for 
extended emission in about ~ 40% of the sources, consistent with the decrease of flux with increasing angular 
resolution. The extended components resemble jet emission in a majority of cases, consistent with the optically 
thin synchrotron emission implied by their steep spectra. We consider the compact 6cm EVN/VLBI radio emission 
of our sources in the context of the "fundamental plane" that previous LLAGN studies identified within the three- 
dimensional parameter space of radio luminosity, X-ray luminosity, and black hole mass. We demonstrate, using 
NGC 7217 and NGC 1068 as particular examples, that high-resolution, multi-epoch radio observations offer useful 
information about the origin of offsets from the fundamental plane. 
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NGC1961 


05:42:04.6 


69:22:43.0 


3934 


49 


SAB(rs)c 


L2 


NGC2782 


09:14:05.1 


40:06:49.0 


2562 


42 


SAB(rs) pec 


H 


NGC3147 


10:16:53.6 


73:24:03.0 


2820 


26 


SA(rs)bc 


S2 


NGC3718 


11:32:35.0 


53:04:04.0 


994 


60 


pec 


L1.9 


NGC4579 


12:37:43.4 


11:49:05.0 


1519 


37 


SAB(rs)b 


S1.9/L1.9 


NGC5953 


15:34:32.4 


15:11:38.0 


1965 


36 


SAa pec 


S2 


NGC7217 


22:07:52.4 


31:21:34.0 


952 


35 


(R)SAB(rs)a 


L2 



Table 1. Basic properties of the NUGA sources taken from NED. a i = inclination of the galaxy. H = HH-region, S = 
Seyfert and L = LINER. 



1. Introduction 

Nuclear activity in galaxies is detected at many different 
levels, supporting the idea that it reflects different phases 
in galaxy evolution, i.e., that every galaxy might exhibit 
nuclear activity at some point in its life. Nuclear activity 
is thought to be due to accretion onto supermassive black 
holes (SMBHs), which are now accepted to reside in vir- 
tually all galaxies, at least those with bulges (Kormendy 
and Richstone 1995; Magorrian et al. 1998). The tight 
correlation of SMBH mass with the stellar velocity dis- 
persion in the bulge (Magorrian et al. 1998; Richstone et 
al. 1998; Gebhardt et al. 2000; Ferrarese & Merritt 2000; 
Tremaine et al. 2002; Shields et al. 2003; Marconi & Hunt 
2003) along with the similar redshift peaks (z w 2) of the 
quasar and starburst populations, suggests a connected, 
or at least similar evolution between starburst and SMBH 
activity. As implied by the population peak, nuclear activ- 
ity seems to increase with redshift, leaving predominantly 
Low-Luminosity AGN (LLAGN) and inactive galaxies in 
the nearby universe. However, many similarities are found 
between LLAGN and High-Luminosity AGN (HLAGN) 
over the entire electromagnetic spectrum, further indicat- 
ing that different activity phases might simply correspond 
to different evolutionary stages. 

Recent studies strongly suggest that HLAGN, 
LLAGN, and even X-ray binaries occupy a so-called "fun- 
damental plane" within the three-dimensional parameter 
space of central black hole mass vs. radio luminosity vs. X- 
ray luminosity (Merloni, Heinz & DiMatteo 2003; Falcke, 
Kording & Markoff 2004). This strongly indicates com- 
mon underlying processes. The exact origin of this cor- 
relation is still the subject of controversy, although it is 
now widely accepted that different accretion rates and/or 
different accretion efficiencies certainly play a significant 
role in setting the various activity levels. In this picture, 
LLAGN may be driven by low accretion rates into a mode 
in which they only release a small fraction of their energy 
through radiation from a radiatively inefficient accretion 
flow (RIAF); the remainder is instead converted into ther- 
mal, turbulent, or mechanical energy. HLAGN, however, 
efficiently cool through radiation, thus producing a high 
activity level. A variety of RIAF models have been invoked 
over the past decade, including Bondi-Hoyle accretion 
(Melia & Falcke 2001), Advection Dominated Accretion 



Flows (ADAFs; e.g., Narayan et al. 1995), Convection 
Dominated Accretion Flows (CDAFs; Ball et al. 2001), 
ADiabatic Inflow Outflow Solutions (ADIOS; Blandford 
& Begelman 1999) and quasi-mono-energetic electron dis- 
tributions (Beckert & Duschl 1997). However, ADAFs and 
Bondi-Hoyle accretion are found to fail for Sgr A* in the 
center of our own galaxy, since they entail gas densities 
which arc significantly higher than observed; they are thus 
already regarded as rather out-of-dat^]]. 

As an alternative to RIAFs, a disk-jet coupling can 
also produce the observed fundamental plane (Heinz & 
Sunyaev 2003); here, the base of the jet rather than a 
RIAF/disk produces the compact (optically thick syn- 
chrotron) emission in LLAGN, while the accretion disk 
dominates the emission in HLAGN (e.g., Falcke & Markoff 
2000). Previous studies of LLAGN (e.g., Ulvestad & 
Wilson 1984; Pedlar et al. 1993; Falcke et al. 2000; 
Nagar et al. 2002, 2005; Anderson, Ulvestad & Ho 2004; 
Anderson & Ulvestad 2005) indicate radio jets and/or 
compact radio emission similar to but less powerful than in 
HLAGN. Jets are intimately connected to accretion onto 
compact objects (e.g., Pringle 1993; Blandford 1993) and 
are thus regarded as a clear signature of accretion-powered 
AGN. However, it is impossible at the moment to distin- 
guish between emission powered by a jet and a RIAF, so 
various combinations of the two must generally be consid- 
ered (e.g., Yuan et al. 2002). In each case, high angular 
resolution observations of the radio emission offer substan- 
tial insight into the nature of the underlying process. 

The seven sources studied in this paper belong to the 
IRAM NU(clei of) GA(laxies) survey (Garcfa-Burillo et 
al. 2003). The NUGA project aims at studying the molec- 
ular gas distributions and kinematics with the highest 
currently available angular resolution and sensitivity in 
a sample of 12 nearby Seyfert and LINER galaxies and 
transition objects (Garcfa-Burillo et al. 2003; Combes et 
al. 2004; Krips et al. 2005; Garcfa-Burillo et al. 2005). 
The study of radio emission is an ideal complement to the 
molecular gas analysis in investigating the activity in these 
sources. We thus selected those sources in our sample with 
existing and clear VLA radio detections (Becker, White & 



1 For a more detailed description of the different RIAFs we 
refer to Quataert et al. (2003), Eckart et al. (2004) and the 
references given above 
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Hclfand 1995; Condon et al. 1996; Ho & Ulvestad 2001) 
for a higher angular resolution study with MERLIN and 
VLBI. 

The paper is organized in the following way. Section [2] 
gives an overview of the observations. Section [3] describes 
the obtained results. The nature of the detected radio 
emission is discussed in Section 2) A comparison of the 
radio luminosities with respect to the fundamental plane 
is presented in Section [5] Finally, the paper ends with a 
summary and conclusions in Section [6] 

2. Observations 

Radio snapshot observations of seven NUGA sources 
(Table QJ have been carried out at 18 cm and 6 cm using 
MERLIN and EVN/VLBI between 2001 and 2004. Two 
sources from this sample (NGC 3718 and NGC 3147) have 
been used for deeper integration (6-9 hours per source). 
The history of the observations is listed in Table [2j We 
selected these seven sources from the NUGA survey be- 
cause they have all been clearly detected on arcsecond 
scales before (VLA observations: Becker, White & Helfand 
1995; Condon et al. 1996; Ho & Ulvestad 2001) but lacked 
consistent mapping at higher angular resolution. With 
both MERLIN and EVN/VLBI, we took advantage of the 
phase-reference technique (Beasley & Conway 1995) by 
observing alternately the target and a phase calibrator 
located within a few degrees. The EVN/VLBI observa- 
tions were conducted with the Mark IV system at a rate 
of 256 Mbs" 1 using 8 IFs, each of 8 MHz bandwidth, 
and 2-bit sampling (mode-256-8-2). The sources were ob- 
served in LCP. At MERLIN, polarization measurements 
were also carried out but without any significant detec- 
tion. The visibility amplitudes were calibrated assuming 
Kelvin to Jansky conversion factors, gain-elevation depen- 
dencies, and system temperatures measured during the 
observations. This and careful amplitude self-calibration 
using the available closure amplitudes provided a typical 
calibration accuracy of 5 — 10%. Approximate beam sizes 
range from 0.05 to 0.5" for MERLIN and 0.001 to 0.04" for 
EVN. RMS image noises range from ~1 /iJy to ~1 mJjd. 
All data were reduced with the AIPS package and mapped 
either in AIPS itself or using DIFMAP (van Moorsel et al. 
1996; Shepherd 1997). 

Additionally, we have information on the continuum 
emission at 3 mm and 1 mm for all seven sources, which 
was gained during the CO(1-0) and CO(2-l) observa- 
tions of the NUGA sample at the IRAM Plateau de Bure 
Interferometer (PdBI). For a more detailed description 
of the PdBI observations, we refer the reader to previ- 
ous NUGA publications (e.g., Garcfa-Burillo et al. 2003; 
Combes et al. 2004; Krips et al. 2005; Garcfa-Burillo et al. 
2005). 

2 snapshot corresponds here to ~3-5 hours on-source inte- 
gration time. 

3 The lower rms noise corresponds to the two deep integra- 
tions, not the snap-shot observations 



3. Results 

We fitted Gaussian profiles in the image plane to all of 
our radio maps to get estimates of component sizes (by 
using the fitted FWHM of the continuum emission and 
half the beamsize as a cutoff lower limit) and to derive the 
respective fluxes ( Table [3] & 01 and positions (Table[3]). We 
first start with some general results for the entire sample 
in Section 13.11 and then discuss each galaxy individually 
in Section I3~2l 

3.1. General results 

Tables [3] and [4] list the derived positions, fluxes and sizes 
of the respective sources. The position of each component 
is consistent across all data sets. Figs. [T] to [3] show the 
continuum maps of the seven sources for each frequency 
and instrument. We find extended emission on all angular 
scales. This is confirmed by comparing the peak intensities 
and integrated flux densities with each other (a difference 
indicating resolved emission). Extended emission is also 
implied by the decrease of flux density as we move from 
VLA to MERLIN to EVN resolution at a fixed frequency 
(e.g., in NGC 1961 and NGC 7217; see Table for the 
VLA luminosities). In at least three objects (NGC 2782, 
NGC 5953, and NGC 7217), we see clear signs of a jet 
on MERLIN scales, and in another two (NGC 1961 and 
NGC 3718) we find indications of jet-like morphologies. 
Interestingly, while for the stronger sources we find flat 
to inverted spectra, the weaker sources have steep spectra 
(Fig. Hand Table El). 

3.2. Results on individual galaxies 

In the next subsections we will present separately the in- 
dividual results for each galaxy. 

We determine for all galaxies the nuclear bolometric 
luminosities via the X-ray flux unless otherwise stated. 
Ulvestad & Ho (2001) propose an empirical relation to es- 
timate the bolometric luminosity of the nucleus in active 
galaxies. Based on a sample of 10 objects, they find Lboi = 
6.7 x L x (2-10keV). Assuming this equation and taking 
the X-ray luminosity (Table 0) published by Roberts & 
Warwick (2000), Ulvestad & Ho (2001), Fabbiano et al. 
(1992), or Terashima et al. (2002), we can crudely esti- 
mate the nuclear bolometric luminosity for each galaxy. 
In some cases, this estimate should be taken as an up- 
per limit: particularly for Roberts & Warwick (2000), a 
significant contribution of circumnuclear star formation 
to the X-ray flux cannot be entirely excluded, due to 
the low (~ 10") angular resolution of the observations. 
Furthermore, we have hard (2-10 keV) X-ray luminosi- 
ties for only three sources, requiring that we extrapolate 
the hard X-ray flux from the soft 0.2-4 keV (or 0.2-2 keV) 
band for the remaining sources as a crude estimate. These 
values are only regarded as rough estimates of the nu- 
clear hard X-ray luminosities, given significant uncertain- 
ties in the contribution of emission from the host galaxy, 
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NGC 5953, NGC 7217 


MERLIN 


4.99GHz 


2002/2003 




all 


EVN 


4.99GHz 


June - 2003 


Eb,Jb,Wb,Mc,Nt, 


NGC 3147, NGC 3718 








On85,Sh,Ur,Tr 




global VLBI 


4.99GHz 6 


May - 2004 


Eb,Jb,Wb,Mc,Nt, 


NGC 1961, NGC 2782, NGC 3718, 



On85,Sh,Ur,Tr,VLBA NGC 5953, NGC 7217 

Table 2. Radio observations for seven NUGA sources. a NGC4579 was excluded since MERLIN observations were 
already carried out; the data are taken from the MERLIN archive. b NGC 3147 and NGC 4579 were not observed 
because VLBA data were already available from Ulvestad & Ho (2001). c The first observations in 2002 were not 
successful due to technical problems with some of the antennas and were thus repeated in 2003; the results of both 
observations are consistent with each other but only the 2003 data sets are used for further discussion due to the 
higher quality of the data. d Eb=Effelsberg, Jb= Jodrell Bank,Wb=Westerbork, Mc=Medicina,Nt=Noto,On85=Onsala 
25m, Sh=Shanghai, Ur=Urumqi,Tr=Torun (see EVN homepage for more information on the respective telescopes: 
http: / / www.evlbi.org/ ) . 



the spectral shape in this band (sometimes a break oc- 
curs at ~1 keV), and the amount of absorption. The last- 
mentioned uncertainty might be the most important one, 
as all of the sources studied in this paper are classified as 
type 2 Seyfert or LINER, implying an obscuration of the 
central engine in the standard AGN unification paradigm. 
However, previous X-ray studies of type 2 Seyfert and 
LINER galaxies, including NGC 3147, NGC 4579, and 
NGC 7217 (e.g., Terashima et al. 2002), indicate that even 
these type 2 LLAGN typically show spectral indices of 
T ~1.6-2.0 (with P(E)=A-E~ r ), suggesting that our ap- 
proximation (r >1) might give a reasonable and reliable 
upper limit. 

We estimate black hole masses for all galaxies using 
the Aibh-Cs correlation, with er s being the stellar veloc- 
ity dispersion of the bulge component (e.g., Gebhardt et 
al. 2000; Tremaine et al. 2002). The velocity dispersion 
has been taken from McElroy et al. (1995) for almost all 
galaxies except for NGC3718, for which we estimated it 
from FWHM(N [II]) taken from Ho, Filippenko & Sargent 
(1997; see also caption of Table [7]). As a caveat, one has 
to keep in mind that the Mbh-c s correlation has been 
mostly established for elliptical galaxies, while our sam- 
ple is based on late type spiral galaxies. However, recent 
results from reverberation mapping indicate that the rela- 
tion still holds for spiral galaxies containing a Seyfert nu- 
cleus (e.g., Wandel et al. 1999; Kaspi et al. 2000; Nelson et 
al. 2004; Onken et al. 2004; Peterson et al. 2005). A differ- 
ent, recently significantly improved approach to determine 
the black hole mass is the bulge-to-black-hole mass corre- 
lation (e.g., Marconi & Hunt 2003, Haring & Rix 2004). 
An estimate of the bulge masses for the NUGA sample is 
in preparation (Hunt et al.) so that we refer the reader 
to this paper, which will also include a comparison of the 
various black hole mass estimates. 



3.2.1. NGC 1961 

NGC 1961 contains a LINER type 2 nucleus in an 
SAB(rs)c host galaxy at a distance of 52 Mpc. This 
is the most distant source in our survey. Assuming 
er s =255km s _1 (McElroy 1995), M bh is inferred to be 
3.OxlO 8 .M (Table [3). Taking L K ~ 4 x 10 40 erg s" 1 
(Roberts & Warwick 2000), the nuclear bolometric lumi- 
nosity of NGC 1961 can be estimated as ^10 41 ergs _1 im- 
plying a sub-EddingtorJf| system with L\, Q \/ L c( \ ~ 10~ 6 — 
10- 5 (Table EJ. 

The MERLIN 18cm map clearly shows extended emis- 
sion, consistent with the difference between peak and in- 
tegrated flux densities (see Fig. [1] and Table [4| and pre- 
vious VLA maps which also indicate extended emission 
(Condon, Anderson & Broderick 1995). The deconvolved 
size is slightly bigger than the beam size, giving further ev- 
idence for extended emission. A two-component Gaussian 
fit assuming a core and a "jet-like" feature results in a 10cr 
determination of the core flux but only in a ~2er estimate 
for the flux of the potential jet. The MERLIN 6cm map 
independently indicates extended emission and, besides 
a compact core, three additional components. The north- 
eastern and the southeastern features lie in the direction of 
the extensions seen in the MERLIN 18cm map. However, 
given that the southwestern feature appears to be sym- 
metrically positioned on the other side of the core with 
respect to the northeastern feature, sidelobe effects can- 
not be totally excluded for these two components. Also, 
since the signal-to-noise ratios of all three "off-core" com- 
ponents are still marginal, they remain questionable. We 
find indications of extended emission in the EVN 18cm 
data as well by comparing the peak with the integrated 
flux density. The distribution appears to be slightly re- 
solved and shows an elongation to the northeast but with 



4 L cd /L Q = 3.22 x 10 4 M/X© 
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Name RA [h:m:s] Aq[s] Dec [°:':"] AS["] 

NGC1961-core 05:42:04.6477 ±0.0003 69:22:42.375 ±0.004 

NGC2782-core 09:14:05.1124 ±0.0003 40:06:49.316 ±0.004 

NGC2782-jet 09:14:05.105 ±0.001 40:06:49.27 ±0.01 

NGC3147-core 10:16:53.6500 ±0.0003 73:24:02.680 ±0.004 

NGC3718-core 11:32:34.8530 ±0.0003 53:04:04.518 ±0.004 

NGC3718-jet 11:32:34.854 ±0.001 53:04:04.523 ±0.004 

NGC4579-core 12:37:43.522 ±0.001 11:49:05.498 ±0.004 

NGC5953-core 15:34:32.383 ±0.002 15:11:37.59 ±0.02 

NGC5953-jet 15:34:32.390 ±0.003 15:11:38.03 ±0.02 

NGC7217-core 22:07:52.3933 ±0.0003 31:21:33.646 ±0.004 

Table 3. Radio positions of the seven NUGA sources, mainly derived from the EVN 6cm (core) or MERLIN 18cm/6cm 
(jet) observations by fitting Gaussian profiles to the emission (the positions are coincident within the uncertainties 
of the different observations). The positional uncertainties are estimated from the Gaussian fits and also include the 
most recent (i.e., determined in the past 4 years) astromectric errors of the used phase-calibrators. 



a different PA than found in the MERLIN maps. The ex- 
tensions are in the direction of the beam, so we must again 
consider sidelobe effects. The two artificial components to 
the west and to the east are found to be at the positions of 
the sidelobes. The 6cm EVN component is only marginally 
resolved and the deconvolved size is smaller than the beam 
size (by a factor of ~2). The PdBI 3mm map also reveals 
a jet-like component of ~lmJy which is ^4" northeast 
of the nucleus (Baker et al. in prep.). Further radio ob- 
servations are needed to clarify the different orientations 
of the jet-like components found at different wavelengths 
and angular resolution. 

Both EVN and MERLIN data imply a spectral index 
of a = —0.3 (for f v cx v a ), indicating a decrease of flux 
with decreasing wavelength. The higher resolutions of the 
6 cm maps mean that the spectral index might be higher 
and so the spectrum even flatter. The upper limit on the 
spectral index derived from the PdBI fluxes is similar to 
the cm-derived values and also suggests a flat spectrum. 

3.2.2. NGC 2782 

The SAB(rs)pec host galaxy of NGC 2782 (D=35 Mpc) 
harbours a dominant nuclear HII region; historically, its 
nuclear activity has been supposed to be mainly domi- 
nated by a powerful nuclear starburst rather than by an 
AGN. However, recent observations of water maser emis- 
sion in this object suggest that a hidden AGN might be 
present as well (Braatz et al. 2004). Its central black 
hole mass can be estimated as 4x10 7 .Mq assuming 
er s =146km s _1 (McElroy 1995). The nuclear bolometric lu- 
minosity amounts to ^10 41 erg s _1 assuming Lx from 
Table [Jj A comparison with the Eddington luminosity 
indicates that NGC 2782 might also radiate at a sub- 
Eddington rate with Lboi/^cd — 10~ 5 — 10~ 4 . 

A strong core and a strong extended jet are visible in 
both the MERLIN 18cm and 6cm maps (see Fig. [1] and 
Table |4|) , supporting the presence of an optically hidden 
AGN; VLA maps also show extended and diffuse emission 
on arcsecond scales (Saikia et al. 1994). The jet is elon- 
gated towards the south/southeast and has a length of 



^0.3" (~ 53 pc). The core seems to be slightly extended, 
too. Going to 6cm and thus to higher angular resolution, 
the core and the jet become more resolved, as indicated 
also by their lower fluxes. The positions of both compo- 
nents agree at the two wavelengths. The EVN 18cm maps 
also show an extended component, again consistent with 
the different peak and integrated flux densities. The ex- 
tension is not in the direction of the beam. However, the 
PA of the extended emission is slightly different from that 
in the MERLIN maps, raising doubts about its reliability. 
In the EVN 6cm map, we find a ~3cr peak at the posi- 
tion of the nucleus. Resolution effects may play a role for 
this tentative detection. The spectral index derived from 
the MERLIN data agrees with the upper limit obtained 
from the EVN data. Parts of the jets might affect the core 
emission in the MERLIN maps since core and jet emission 
cannot be easily separated from each other in the 18cm 
MERLIN map, complicating the estimate of the spectral 
index for the core and the jet. 

No continuum emission was detected at 3 mm or 1 mm 
with the IRAM PdBI. The derived upper limits are given 
in Tabled] 

3.2.3. NGC 3147 

NGC 3147 is classified as SA(rs)bc galaxy at a distance 
of 38 Mpc with a Seyfert type 2 nucleus. The central 
black hole mass can be derived as ^4xl0 8 .M© with 
er s = 268km s _1 (McElroy et al. 1995) and the nuclear 
bolometric luminosity is ^2xl0 42 erg s _1 (Table [7]), im- 
plying L hol /L ed ~ 5x 10" 5 . Therefore, NGC 3147 radiates 
at a sub-Eddington rate. 

NGC 3147 is among the strongest radio sources in our 
survey (see Fig.Q]and Tabled]). In all maps (including the 
mm-data and VLA data (Ho & Ulvestad 2001), it appears 
to be point-like. The almost identical fluxes determined 
at both wavelengths and instruments are a further sign of 
very compact emission in this object. The flat to inverted 
spectral index derived from the EVN/MERLIN flux den- 
sities indicates that no extended jet is present. The point- 
like structure and the spectral indices (at cm- wavelengths) 
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Table 4. Results of the radio and mm continuum observations with MERLIN, EVN and PdBI. a from the MERLIN 
archive (PI: N.Nagar). b taken from VLBA observation carried out by Ulvestad & Ho (2001; =Ulv01). Please note that 
the flux errors do not include the calibration uncertainties of 5-10% and represent only the statistical uncertainties of 
the fit. The given angles are derived from North to East. 



are consistent with what Ulvestad & Ho (2002) found 
in VLBA data (see also Anderson, Ulvestad & Ho 2004; 
Anderson & Ulvestad 2005). Our two observing epochs for 
the EVN 18cm data show no evidence for any variability 
of this source. We find a turnover between the cm- and 
the mm- fluxes (Fig. B} , as can also be seen in the differ- 
ent spectral indices at mm- and cm- wavelengths (Table [5]). 
The turnover is most likely caused by synchrotron self ab- 
sorption, as is discussed in further detail in a separate 
paper (Krips et al. 2006). 



3.2.4. NGC 3718 

The galaxy NGC 3718 is located at a distance of 13 Mpc, 
hosts a LINER type 1.9 nucleus, and has a warped gas 
and dust disk which might be caused by an interaction 
with a close companion galaxy (e.g., Krips et al. 2005). 
A black hole mass of ^5xlO 7 .M0 is implied by a s = 
157km s _1 (Ho, Filippenko, & Sargent 1997), and a nu- 
clear bolometric luminosity of ~10 41 erg s~ x (Table [7]) im- 
plies a sub-Eddington system with Lboi/^ed — 10 -5 — 
10~ 4 . 

Similar to NGC 3147, NGC 3718 is among the 
strongest radio continuum sources in our survey (Table |4]). 
It shows extended emission on MERLIN 18cm scales and 
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Fig. 1. MERLIN (2 left columns) and EVN/VLBI (two right columns) of the 18cm and 6cm continuum emission in 
NGC1961 (upper panels), NGC2782 (middle panels) and NGC3147 (lower panels). The respective synthesized beams 
are shown in the lower left part of each map. Contour levels for MERLIN-18cm (NGC 1961; upper left): (3<r=)0.57 
to 2.28mJy/beam in steps of la; for MERLIN-6cm (upper left middle): (3<7=)0.39 to 0.78mJy/beam in steps of la; 
EVN-18cm (upper right middle): (3cr=)0.21 to 0.49mJy/beam in steps of la; EVN-6cm (upper right): (3cr=)0.14 to 
0.61mJy/beam in steps of lcr; MERLIN-18cm (NGC2782; middle left): (3er=)0.45 to 1.7mJy/beam in steps of lcr, 
for MERLIN-6cm (middle left middle): (3cx=)0.39 to 0.78mJy/beam in steps of Ict, for EVN-18cm (middle right 
middle): (3o=)0.18 to 0.36mJy/beam in steps of lcr; for EVN-6cm (middle left): (2<t=)0.2 to 0.3mJy/beam in steps 
of la. MERLIN- 18cm (NGC3147; lower left): (5er=)0.75 to 7.5mJy/beam in steps of 5cr, for MERLIN-6cm (lower left 
middle): (3o=)0.72 to lO.lmJy/beam in steps of 5a, for EVN-18cm (lower right middle): (3o=)0.65 to 5.2mJy/beam 
in steps of 5cr; for EVN-6cm (lower left): (5o=)1.4 to 8.3mJy/beam in steps of 5a. Negative contours, if visible, 
correspond to 3<r. 



signs of a quite compact (~0.5" =34 pc) 4cr jet in a north- 
west direction (Fig. [2]). This is supported by VLA maps 
that also reveal extended emission on large scales (Condon 
1987). There are some weak indications for extended emis- 
sion from inspection of the 6cm maps and comparison of 
the peak and integrated flux densities. Moreover, the EVN 
6cm map reveals a slight extension but, in contrast to the 



MERLIN 18cm map, in a westward direction, which we 
therefore do not identify as the same jet. The flux den- 
sities from MERLIN and EVN are identical, indicating 
that there is no emission extended on intermediate scales 
in NGC 3718 other than the small-scale jet. Our fluxes are 
in good agreement with those obtained in VLBA obser- 
vations by Nagar et al. (2002). The spectral indices are 
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Fig. 2. MERLIN (2 images on the left respectively) and EVN (2 images on the right respectively) of the 18cm and 
6cm continuum emission in NGC3718 (upper panel), NGC4579 (middle panel) and NGC5953 (lower panel) . The 
respective synthesized beams are shown in the lower left corner of each map. MERLIN-18cm (NGC 3718; upper 
left): (5er=)0.64 to 5.1mJy/beam in steps of 5a, for MERLIN-6cm (upper left middle): (5a=)0.5 to 5.0mJy/beam 
in steps of 5a, for EVN-18cm (upper right middle): (5a=)0.75 to 4.5mJy/beam in steps of 5a; for EVN-6cm (upper 
right): (5<r=)0.55 to 5.0mJy/beam in steps of 5a; Contour levels for MERLIN-6cm (NGC4579; middle left middle): 
(5cr=)1.0 to 15.0mJy/beam in steps of 10a, for EVN-18cm (middle right middle): (3a=)3.3 to 16.5mJy/beam in steps 
of 3ct; MERLIN- 18cm (NGC5953; lower left): (3a=)0.48 to 1.5mJy/beam in steps of la, for MERLIN- 6cm (lower left 
middle): (3a=)0.74 to 1.7mJy/beam in steps of la, for EVN-18cm (lower right middle): (2a=)0.12 to 0.24mJy/beam 
in steps of la; for EVN-6cm (lower right): (2a=)0.2 to 0.4mJy/beam in steps of la. Negative contours, if visible, 
correspond to 3a. 



very consistent for the MERLIN and EVN data and are 
derived to be 0.4, describing an inverted spectrum. 

3.2.5. NGC 4579 

NGC 4579 is a SAB(rs)b galaxy at a distance of 
20 Mpc with a Seyfert type 1.9 and/or LINER type 



1.9 nucleus. The black hole mass is ~10 8 Mq assuming 
a s =185km s-^McElroy 1995). Ho et al. (1999) determine 
the nuclear bolometric luminosity to be ~10 42 erg s -1 . 
This gives a value of ~10 -4 for Lboi/^od indicating a sub- 
Eddington system. 

NGC 4579 was already observed with MERLIN at 
18cm (PI: N.Nagar, from the MERLIN archive) and with 
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Fig. 3. MERLIN (2 images on the left) and EVN (2 images on the right) of the 18cm and 6cm continuum emission 
in NGC7217. The respective synthesized beams are shown in the lower middle part of each map. Contour levels 
for MERLIN-18cm (left): (3cr=)0.36 to 1.4mJy/beam in steps of la, for MERLIN-6cm (left middle): (3cr=)1.5 to 
4.0mJy/beam in steps of la, for EVN-18cm (right middle): (3cr=)0.21 to 0.5mJy/beam in steps of la; for EVN-6cm 
(right): (3<r=)0.2 to 0.6mJy/beam in steps of la. Negative contours, if visible, correspond to 3a. 



the VLBA at 18cm, 13cm, 6cm and 4cm (Ulvestad & Ho 
2001). We thus conducted only complementary observa- 
tions with MERLIN at 6cm and EVN at 18cm (Fig[2]and 
Tabled]). In the VLBA maps from Ulvestad & Ho (2001), 
no extended emission can be seen, while VLA maps from 
Ho & Ulvestad (2001) at 20cm and 6cm show a jet-like 
extension of —3" to the northwest. Our 18 cm and 6 cm 
maps do not show any jet-like feature on sub-arcsecond 
scales, but the difference between peak flux and integrated 
flux still indicates extended emission. Furthermore, the de- 
convolved size at 18cm is slightly larger than the beam, 
supporting the presence of extended emission. Falcke et 
al. (2001) report variability of the radio flux at 15 GHz on 
timescales of 1-3 years. We too find a difference of ~ 10% 
between the EVN and VLBA 18 cm flux densities. Thus, 
the determination of the spectral index must be made and 
interpreted with caution since not all fluxes were measured 
simultaneously. This might explain the large discrepancy 
between the MERLIN 18cm flux and that from the EVN 
18cm data, which is higher by a factor of 2-3. 

3.2.6. NGC 5953 

NGC 5953 is characterised as a SAa pec galaxy at a dis- 
tance of 26 Mpc hosting a Seyfert type 2 nucleus. It is also 
known to interact with a close companion NGC 5954 (e.g., 
Iono, Yu & Ho 2005). The black hole mass of ~7x 10 6 A^ Q 
(cr s =94km s _1 ; McElroy 1995) is the smallest one m our 
sample. Woo & Urry (2002) find a surprisingly high nu- 
clear bolometric luminosity of — 10 44 erg s _1 but their es- 
timate is based on fluxes taken over very large apertures. 
No X-ray luminosity has been published yet. Thus, any 
conclusion about its accretion efficiency remains specula- 
tive. 

NGC 5953 turns out to be the weakest radio source in 
our sample (Tabled]). It is clearly detected in the MERLIN 
18cm and 6 cm maps, while it remains undetected in the 
EVN maps (Fig. [3]) . We find 4cr peaks in the latter at the 
position of the nucleus. The MERLIN 18cm data show a 
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Table 5. Spectral indices of the respective sources. S v oc 



taken from Ulvestad & Ho 2001. 
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NGC1961 


5.0 
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>4xl0 Y 


NGC2782 
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Table 6. Brightness temperatures estimated with the for- 
mula given in Condon et al. (1982). The frequency (v), the 
flux density (S v ) and the deconvolved size (9) as an upper 
limit are taken from Tabled] a taken from Ulvestad & Ho 
2001. 



strong jet to the north, which is slightly resolved with the 
MERLIN beam (—0.3" =41pc). A steep spectral index for 
the core is derived from the MERLIN data, which might 
be biased by resolution effects in the MERLIN 6cm map. 
VLA maps (e.g., Condon et al. 1990; Iono, Yu & Ho 2005) 
also indicate extended emission on large angular scales, 
but the emission is somewhat smeared between the two 
interacting systems in the VLA beam. 
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Fig. 4. Spectral energy distribution from the 
MERLIN/PdBI {upper panel) and EVN data {lower 
panel). The dashed line in the upper panel separates the 
MERLIN from the PdBI data. Values for SgrA* were 
taken from Falcke et al. 1998 and Zhao et al. 2001. 



3.2.7. NGC 7217 

NGC 7217 contains a LINER type 2 nucleus in an 
(R)SAB(rs)a host galaxy at a distance of 13 Mpc. It 
has a black hole mass of ~3xl0 7 M.® assuming a s = 
132km s (McElroy 1995) and an estimated nuclear bolo- 
metric luminosity of ~10 40 erg s _1 (Table [7J). Both values 
result in L bol /L ed ^ 10~ 6 - 10~ 5 . Thus, NGC 7217 is also 
a sub-Eddington system. 

The MERLIN 18cm map (Fig. [3] and Table SJ clearly 
shows an extended component with an elongation to the 
south, although NGC 7217 appears to be unresolved on 
VLA scales. The MERLIN 6cm and EVN 18cm compo- 
nent remains pointlikc, however, while the EVN 6cm map 
reveals an extended component but in an eastern direc- 
tion different from that in the MERLIN 18cm map. This 
apparent inconsistency might be a result of the different 
sensitivity levels and angular resolutions obtained by the 
respective observations. The MERLIN and EVN spectra 
are the most strongly inverted in our sample. 



4. Nature of the radio emission: thermal or 
non-thermal? 

Based on the equation given by Condon et al. (1982) and 
taking the deconvolved source size as an upper limit for the 
actual source size and the flux density at 6 cm from Table|4] 
(except for NGC 2782), we can estimate lower limits for 
the brightness temperature in each of the seven galaxies. 
For most of them, the brightness temperatures are high 
enough (T;, > 10 s K) to be consistent with synchrotron 
emission given that the assumed source sizes are only up- 
per limits. For NGC 3147, for instance, Anderson et al. 
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Fig. 5. Comparison of radio (5 GHz) with X-ray luminosi- 
ties {upper panel) and with black hole mass {lower panel:). 
We extrapolated the 15 GHz literature data {triangles) to 

5 GHz where necessary using the spectral indices indi- 
cated in Merloni, Heinz & Di Matteo (2003); this might 
cause a larger scatter of the data because of their larger 
uncertainty. The filled light grey triangles represent VLA 
data while the filled dark grey triangles are taken with the 
VLBA. For those galaxies from our sample, which lack 
published hard (2-10keV) X-ray luminosities {open black 
circles), we assumed Lx(0.2-4keV)>Lx(2-10keV) as sim- 
plification (see text for discussion). The open grey circle 
represents NGC 5953. The literature values are taken from 
Terashima et al. (2002; X-ray), Nagar, Falcke & Wilson 
(2005; VLBA or VLA data at 15 GHz) and Merloni, Heinz 

6 Di Matteo et al. (2003; M BH )- Errors of the radio lu- 
minosities are from the EVN 5.0 GHz data. The radio 
data {blue star) from NGC 1068 refer to component SI 
which is supposed to be the AGN and are taken from 
Gallimore et al. (2004), while the black hole mass is from 
Hure (2002) and Lodato & Bertin (2003). The X-ray lumi- 
nosity is adopted from Merloni, Heinz & Di Matteo (2003) 
but has to be regarded as a lower limit due to a possible 
heavy obscuration of the AGN in NGC 1068. 
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Table 7. Black hole masses, bolometric and Eddington luminosities for the seven galaxies. a assuming Mbh = (1.3 ± 
0.1) x IOIMqX (<j s /200[km s -i]) 4 0±0 - 3 f rom Tremaine et al. 2002 with <r s taken from McElroy et al. 1995. b a s taken 
from Ho, Filippcnko & Sargent (1997) and assuming «7 S =FWHM([0 III])/V8 In 2=FWHM( [N 11]) /V& In 2. c taken 
from Roberts & Warwick 2000. d taken from Ulvestad & Ho 2001. e taken from Fabbiano et al. 1992. ' taken from 
Terashima et al. 2002. 9 assuming Lb o i=6.7xi x (2-10keV) from Ulvestad & Ho 2001; as simplification and since we were 
interested only in a rough estimate of the bolometric luminosity, we assumed L x (2-10keV)« L x (0.2-4keV) allowing for 
half an order of magnitude uncertainty. h taken from Woo & Urry (2002). 1 taking the standard Eddington equation. 
J assuming L5.0 = ^5.0 x 47rD 2 S'„ where S„ is the integrated flux density from Table HI k Ho et al. 1999. 1 derived 
via D~V hc i/74Mpc. m extrapolated from 15 GHz VLA observations by Nagar et al. (2005) except for NGC 2782, 
NGC 5953 and NGC 7217 which were taken from 1.4 GHz VLA observations done by Becker, White & Helfand (1995), 
Condon et al. (1996) or Iono, Yun & Ho (2005). n We corrected the X-ray luminosities to match the distances used in 
this paper, if significantly different luminosity distances were assumed in the X-ray data papers. 
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Fig. 6. Comparison of radio fluxes with a projection of 
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Di Matteo (2003). See caption of Fig. [5] for more details 
on data. 



(2004) derive a source size of <1 mas, increasing the tem- 
perature determined from our data to ~4xl0 8 K. Given 
the loose upper limit on the source size for NGC 2782, its 
temperature might also be compatible with synchrotron 
emission. However, the lower limit on virtually ex- 
cludes the possibility that the radio emission in the cen- 
tral ±10 mas of NGC 2782 is coming from star formation, 
which typically produces < 10 5 K, supporting the presence 



of synchrotron (self-absorbed) emission from an optically 
hidden AGN. 

The situation for NGC 5953 and NGC 7217 might 
be different. Although NGC 5953 has a very high up- 
per limit on the source size, its brightness temperature 
would only go up to ~10 7 K when taking a source size of 
1 mas, probably still too low to agree with synchrotron 
self- absorption from a compact source but also too high 
for star formation. However, optically thin synchrotron 
emission from an (extended) jet might not disagree with 
the low brightness temperature; this would also fit to the 
steep radio spectrum found in NGC 5953. The continuum 
emission in NGC 7217 appears already slightly resolved 
in the VLBI 6 cm map, leaving little leverage to increase 
the true brightness temperature to 10 s K excluding thus 
a self-absorbed compact synchrotron source. In contrast 
to NGC5953, the inverted radio spectra of NGC7217 are 
not in agreement with optically thin synchrotron emis- 
sion. Alternatively, thermal free-free emission or electron 
scattered synchrotron emission may be the mechanism re- 
sponsible for the continuum similar to NGC 1068 (e.g., 
Roy et al. 1998; Gallimore et al. 1997, 2004; Krips et al. 
2006). Although we see no evidence that starbursts dom- 
inate the centers of our LLAGN, the larger VLA fluxes 
imply the presence of larger scale diffuse emission that 
might be produced by starbursts. 

5. The fundamental plane 

A highly significant relationship linking X-ray luminos- 
ity, black hole mass, and radio luminosity has been found 
to hold for many black hole systems, including both X- 
ray binaries and active galaxies (e.g., Merloni, Heinz & 
Di Matteo 2003; Falcke, Kording & Markoff 2004). In 
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Fig. [5] & El the VLBI luminosities at 5 GHz (L5.0) are 
compared with published X-ray luminosities (Fabbiano et 
al. 1992; Roberts & Warwick 2000; Terashima et al. 2002) 
and black hole masses (McElroy et al. 1995) for six of 
our objects, along with literature data for —30 LINER 
and Seyfert galaxies taken from Terashima et al. (2002), 
Merloni, Heinz & Di Matteo (2003) and Nagar, Falcke & 
Wilson (2005). In contrast to Merloni, Heinz & Di Matteo 
(2003), and Falcke, Kording & Markoff 2004, -50% of 
the plotted radio fluxes were obtained with the VLBA or 
VLBI, i.e., at higher angular resolution, reducing the con- 
tamination of nuclear emission by extended components 
such as (large-scale) jets or star formation. 

We find that our new data are overall in good agree- 
ment with the existing data and lie within the same em- 
pirical "fundamental plane" derived by Merloni, Heinz & 
Di Matteo (2003) and Falcke, Kording & Markoff (2004), 
even in those cases for which we used a crude estimate for 
the hard X-ray luminosity. This relation has been inter- 
preted theoretically in terms of a disk-jet coupling, as sug- 
gested by Heinz & Sunyaev (2003), radiatively inefficient 
accretion flows, or a combination of both. Radiatively 
efficient accretion flows can be almost completely ruled 
out for these sub-Eddington systems (Merloni, Heinz & 
Di Matteo 2003). The larger scatter in the literature 
data compared to Heinz & Di Matteo (2003) and Falcke, 
Kording & Markoff (2004) might be caused by the uncer- 
tainty of the assumed spectral indices that were used to 
extrapolate the 5 GHz fluxes from the published 15 GHz 
data. 

Within our sample, NGC 7217 seems to have the 
largest offset from the fundamental plane. As mentioned 
in Section [4] NGC 7217 does not appear to be consis- 
tent with (direct) synchrotron emission. It is unclear at 
this point whether its emission is produced by electron- 
scattered synchrotron emission, as might be the case in 
NGC 1068 (e.g., Roy et al. 1998; Gallimore et al. 1997, 
2004; Krips et al. 2006), or by thermal free-free radia- 
tion. The large uncertainty of the X-ray luminosity in 
NGC 1068 due to the heavy obscuration of its AGN and 
its large variablity (see discussion in Gallimore et al. 2004) 
also questions the compatibility of NGC 1068 with the 
fundamental plane (Fig. & 0, when the VLBA flux of 
the AGN (component SI) is plotted instead of the lower 
angular resolution VLA flux adopted by Merloni, Heinz & 
Di Matteo (2003). The latter is certainly contaminated by 
emission from the extended parts of the jet. 

6. Summary & Discussion 

This paper presents the first results of snapshot radio ob- 
servations with MERLIN and EVN of a sample of seven 
LLAGN taken from the NUGA survey. Besides compact 
emission, small-scale extended emission in the form of 
jets is also found in these weak radio sources. The radio 
fluxes decrease significantly with increasing angular res- 
olution at the respective wavelength, indicating the pres- 
ence of extended and diffuse emission. With one exception 



(NGC 5953), the radio core components have flat to in- 
verted spectra within the errors. A comparison of the de- 
rived radio luminosities to the black hole mass, estimated 
via the respective stellar velocity dispersions, and the X- 
ray luminosities suggests correlations among these quan- 
tities even at high angular resolution, thus supporting the 
existence of a fundamental plane. However, as NGC 1068 
and NGC 7217 demonstrate, this fundamental plane could 
be violated if the radio emission is produced not by direct 
synchrotron processes, but rather by electron scattered 
synchrotron or free-free emission, even though the latter 
are also expected to be closely connected to X-ray activity. 
Whether or not such discrepancies reflect substantially dif- 
ferent activity mechanisms in NGC 1068, NGC 7217, and 
the remaining LLAGN in our sample has to be further 
investigated. This result emphasizes the need for high an- 
gular resolution observations of the radio emission in these 
sources to verify the validity of the fundamental plane and 
test its limits. 
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